Scaling theory of excitations in trapped fermionic condensates 
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• Formulas are derived for the coupled quadrupolar and monopolar oscillations of a fermion conden- 

' sate trapped in a axially symmetric harmonic potential. We consider two-component condensates 

, with a large particle-particle scattering length and in the limits of frequencies either small (adia- 

, ^ . batic) or large (diabatic) compared to the pairing gap. Using the method of coordinate scaling, 

' we reproduce Stringari's results for adiabatic excitations and present new results for the diabatic 

^ limit. Unlike the adiabatic oscillations, the diabatic ones show sensitivity to the interaction between 

. particles. 
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I. INTRODUCTION 



^ i' Experimental techniques in atomic cooling and trapping have now evolved to a point where condensates of fermionic 
' [ atoms such as ^Li can be produced and studied 0. One of the interesting measurable properties is the spectrum of 
^ ' collective excitations. Stringari 2] has investigated this question for the case of the unitary limit of a two-component 
O ] system. By "unitary" is meant that the scattering length associated with the two-body interaction is large compared 
to all other scales in system. In ref. ^ it is assumed that the system behaves hydrodynamically. This requires the 
dynamics to maintain an isotropic pressure tensor when the system is deformed. The condition for this adiabatic limit 
^ ,: to apply to fermion systems is that they be superfluid and that the collective oscillation frequency be much less than 
^ ■ the pairing gap. The other extreme, the diabatic limit, obtains if the oscillation frequency is much greater than the 
pairing gap. This in fact is the case for the giant resonances in nuclear spectroscopy. In that limit, the dynamics is 
f~i ■ governed by the mean field acting on the particle orbitals and coupling through the ordinary density. We believe it is 
Qh' of interest to investigate this limit for the atomic condensates as well, to understand the dependence of the frequencies 
on all possible parameters of the system. 
I The time-dependent mean-field theory may solved in the small-amplitude limit by the RPA equations, but inclusion 
J> of all the single-particle degrees of freedoms requires a considerable numerical effort. This is especially the case for the 
lO theory including the pairing degrees of freedom as well 3, 4J. Much simpler analytic formulas can be derived under 
, the assumption that the excitations are fully collective, meaning that they exhaust the sum rule for some one-body 
field. In the sum rule approach, the formulas express the oscillation frequency in terms of a ratio of two integrals 
involving only ground-state quantities in the integrands. Again referring to the experience in nuclear physics, this 
approach has been very successful, particularly with the monopole and quadrupole modes of oscillation. 
I Considering atomic condensates in an axially symmetric harmonic trap, it is natural to consider compressional and 
. quadrupolar modes generated by the fields Q = x^,y'^, and z^. These fields induce scaling transformations on the 
^ ■ orbital wave functions, e.g., for the field a;^, the transformed wave function is |^ 

0(x,y,z) ^ = e-l^^/^e-^^x^y^z) (1) 

^' . . . 

I where Px is a scaling coordinate. The collective frequencies may be calculated by constructing the energy function in 

^ ' the scaling coordinates and solving the classical equation of motion. This route is rather transparent and we follow it 
• ^ in the derivation below. Consider an energy function of the form 

H ; E^T + Vext + U (2) 

_ _ I 

where T is the kinetic energy operator, U is the interaction energy, and Vext is the trap field. We write the trap field 
as 

V,.t = l{x' + y' + {\zf), (3) 

setting the mass of the particles and the transverse harmonic oscillator frequency to one. In the Thomas-Fcrmi 
approximation of a free fermion gas the ground state density is given by 

3/2 



Po(r) = c(l-(a;/i?)2-(y/i?)2_(Az/i?)2) , (4) 
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where c is a constant and R is the transverse radius of the condensate. The kinetic energy is given by 

{T)tf^ Jd^pl/'^iSn'f (5) 

in a two-component Fermi gas. The density also has the form of eq. (0J in an interacting two-component fermion 
system in the unitary hmit. The scaled density is given by 

p^(f) = ce-^^-l^y-'^^ (l - {e-^'^x/Rf - {e-l^yy/Rf - {Xer^' z / Rf'^'^'^ . (6) 

We next evaluate the various terms in the energy as a function of the scaling variables. The expectation of the external 
potential is given by 

where the prefactor a is the integral a = J (Prp^x^ /2. The kinetic energy may vary in two different ways, depending 
on whether or not the Fermi surface remains spherical under the dynamic evolution. In the adiabatic limit we have 

{T)ad - 6e-2(/3^+/5«+/3^)/3 (8) 

where b — {T)tf- In the diabatic limit the contributions in the three Cartesian directions are independent and the 
kinetic energy is 

{T),.a^l{e-'^'^+e''f'y+e~'^^). (9) 

We next turn to the potential energy. It behaves very simply in the unitary limit, being proportional to the adiabatic 
kinetic energy with a universal constant of proportionality. We shall further assume that the potential energy behave 
as a scalar, creating an isotropic pressure. It thus depends only on {T)ad- We parameterize it as 

{U)^{i^\){T)ad (10) 

We may now calculate the vibrations in a harmonic approximation making a second-order Taylor expansion of the 
energy function of the collective variables. The restoring force constants are given by the hessian matrix K, with 

We will also need the inertia tensor I associated the collective variables. Defining the velocity field v = /3VQ/2, the 
element 7^ is the coefficient of (3i(3j in the expression J (Prpov ■ v. In our case the integral is 



lij = d rportrj. (12) 



Note that the integral is essentially identical to that for (Vext)- Thus 



1 

I = 2a ( 1 I . (13) 




Ku = uj^lu (14) 

Before applying the formulas for the diabatic frequencies, we confirm previous results on the adiabatic frequencies 0, 
1^. Here we use eq. ^ for the kinetic energy and eq. H10|l for the interaction energy, i.e. setting (E) — £,{T)ad in 
eq. (11). We also use the fact that /3 = in equilibrium to eliminate the coefficient b in the expression for the kinetic 
energy. The equilibrium condition 



(15) 
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is satisfied by setting b = 3a/^. The squared frequencies may tlien be obtained from the eigenvalues of the matrix 



I-V2k,,I-V2 = ^ I 1 4 A I . (16) 



4 1 A 
1 4 A 
A A 4A2 

Note that aU dependence on the interaction strength ^ has disappeared. The eigenvalues of eq. Hlt)|) are 



UJ 



I2 = (5 + 4A2 ± v/25 - 32A2 + 16A4)/3 (17) 



Lol = 2. 

This agrees ref. as well as the frequencies calculated for noninteracting fermions in ref. |^ . To get more insight into 
the physics, let us specialize to the case of a spherical trap (A = 1). The frequencies become ujf = 4, corresponding 
to the monopole vibration and ^1 3 = 2, corresponding to two of the five degenerate quadrupolar vibrations^. 

Next we carry out the same calculation using the diabatic energy. The self-consistency condition eq. (|15|l together 
with the strong-coupling relation eq. determines the internal energy function to be 

{T)d^a + (U) = |(e-2^- + e--'^y + e-2A) + t^^_2l^-2iP.+f^y+M/3 ^ (18) 

Evaluation of the derivatives for the force matrix then yields 

^d^a = ^ f -1 2 -1 I + ^ I 1 4 1 I . (19) 

The eigenvalues of I^-^Z-^K^ial^^^^ are 

LU 




= (/ ± V/2-72eA2(l + 0)/3e (20) 



2 ^ 2 

^3=2+7 

where f = 1 + 2X^ + 5^-1- 'iX^S,- For the convenience of the reader, a script is given in |^ of the mathematical 
operations leading to eqs. H17II20() . 

The frequencies in eq. H2U|I depend on the interaction constant ^, in contrast to the adiabatic frequencies. Their 
measurement would thus provide information about the adiabaticity and thus about the relative size of the pairing 
gap and the vibrational frequency. As a final remark, let us specialize again to the spherical trap. The frequencies 
become uif = A and a;| 3 = 2 + 2/^. The monopole frequency coi is the same as in the adiabatic formula, which 
is not surprising because the Fermi surface remains spherical under the compression by the monopole scaling field. 
The quadrupole modes depend on the interaction and lie above the monopole at the physical value of ^, calculated 
to be ^ « 0.44 with pairing and ^ w 0.55 without pairing 9]. This is completely at variance which the behavior of 
the quadrupole mode in other systems such as atomic nuclei, where the incompressibility of the medium produces a 
higher monopole frequency. 



It is amusing to note that the diabatic harmonic formula for the the giant nuclear quadrupole resonance is also reduced from the nominal 
value of LJ'^ = 4 by a factor of two0,|3|. In both cases, the flow is incompressible so that the co- moving density remains the same during 
the oscillation. In the nuclear case with a short-range interaction, this implies a vanishing interaction contribution to the restoring force 
constant, leaving a kinetic contribution. The situation is converse for the trapped fermion case: here the kinetic contribution vanishes 
leaving only the interaction with the external field. 
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